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In recent years there has been a strong interest in non-perturbative approaches to continuum Yang-Mills theory and 
much research has been carried out in this field. Most of these approaches rely on gauge fixing by the Faddeev-Popov 
method. This concerns, in particular, the Dyson-Schwinger equation approach in Landau gauge (for recent reviews 
see [1, 2]) and in Coulomb gauge [3, 4], and the Hamiltonian approach in Coulomb gauge [5, 6]. Furthermore, in 
Coulomb and Landau gauge the Faddeev-Popov determinant is assumed to dominate the infrared sector of the theory, 
see e.g. Ref. [7]. 

The standard way to renormalize the theory is to represent the Faddeev-Popov determinant by ghost fields and 
use the BRST invariance of the resulting local field theory. In some cases it is, however, advantageous not to 
introduce ghost fields but to keep the Faddeev-Popov determinant explicit and treat it non-perturbatively. This is, in 
particular, the case in the Hamiltonian approach in Coulomb gauge [6]. The renormalization of this theory requires 
then to identify the ultraviolet singular pieces of the Faddeev-Popov determinant and to remove them by appropriate 
counterterms. In Ref. [8] the UV counterterms of the Faddeev-Popov determinant in d = 3 + 1 Coulomb gauge were 
identified using an approximate expression for the Faddeev-Popov determinant [9] 



where \ is the so-called curvature, which represents the ghost loop contribution to the gluon self-energy. The repre- 
sentation (1) can be shown to be correct to 2-loop order for the energy. Not surprisingly the required counterterm to 
lnDet(— Dd) was found from this representation to be given by a mass-like term 



with the coefficient C being linear in the UV-cutoff. In the present paper we will calculate the counterterm to 
the Faddeev-Popov determinant exactly without resorting to the representation (1) and also determine the precise 
numerical coefficient of the counterterm. For this purpose we carry out a heat-kernel expansion of the Faddeev-Popov 
determinant in d = 3,4 Landau gauge. The d = 3 Landau gauge corresponds to the d = 3 + 1 Coulomb gauge, the 
case of most interest to us. 

The organization of the paper is as follows: In Sec. II we briefly summarize the heat-kernel expansion of functional 
determinants. In Sec. Ill we carry out the heat-kernel expansion for the Fadeev-Popov operator in Landau gauge. The 
corresponding heat coefficients are evaluated in Sec. IV and the UV-divergent terms are isolated in Sec. V. Finally in 
Sec. VI we present our conclusions. 



We essentially follow Ref. [10] here. The determinant of a positive definite operator M can be represented by the 
following proper-time integral 



I. INTRODUCTION 
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II. HEAT KERNEL EXPANSION OF FUNCTIONAL DETERMINANTS 




(3) 
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where A is an ultraviolet momentum cutoff and 

K (r) = e' MT (4) 

is the heat-kernel. Furthermore, Tr includes both the trace over the discrete indices tr and the integration over space 
time / d d x. The heat-kernel (4) satisfies the heat equation 

d T K(r) + MK(t) = (5) 

and the boundary condition 

K(t = 0) = 1 . (6) 

In the usual situation the operator under interest M contains a "free" part Mo, which can be treated exactly. In this 
case it is convenient to express the full heat-kernel as 

(x\K(r)\y) = (x\K (r)\y) (x\H(r)\y) , (7) 

where 

K (t) = e~™° (8) 

is the free heat-kernel satisfying 

(d T + M ) K (t) = , K (t = 0) = 1, (9) 

and H(t) embodies all the "interactions". From the boundary conditions to the heat-kernels follows the boundary 
condition for the interaction part 

H(t = 0) = 1. (10) 
The interaction part H(t) is then expanded in powers of the proper time r 

oo 

(x\H(r)\y)=J2hn(x,y)r n , (11) 

where the h n (x,y) are referred to as heat coefficients. Eq. (10) implies 

h (x,x) = l. (12) 
With (7) and (11) we find for the functional determinant (3) 



n-1 



TrlnM = -V/ dr r™" 1 Tr (K (T)h n ) . (13) 



The advantage of the heat-kernel expansion (13) is that with increasing order n the terms become less and less 
ultraviolet singular and only the first few terms are ultraviolet singular while, starting at a certain n (which depends 
on the number of space-time dimensions), the terms are ultraviolet finite. In the following we shall apply this heat- 
kernel expansion to the Faddeev-Popov determinant in Landau gauge in d = 3, 4 Euclidean dimensions. 

III. HEAT KERNEL EVALUATION OF THE FADDEEV-POPOV DETERMINANT IN LANDAU 

GAUGE 

The Faddev-Popov kernel in Landau gauge 

dA = (14) 

is given by 

M — —Dd , (15) 
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where 



D = d + A, A = A a T a , (T a 



,6c 



bac 



(16) 



denotes the covariant derivative in the adjoint representation of the gauge field. Here and in the following all matrix 
valued quantities will be defined in the adjoint representation of the gauge group. Furthermore, we have absorbed 
the coupling constant in the gauge field. Note that the Faddeev-Popov operator (15) is positive definite in the first 
Gribov region, to which the gauge fields should be restricted, and hence its determinant can be expressed by the 
proper-time integral (3). One should mention, however, that the Faddeev-Popov operator has constant zero modes 
corresponding to the global gauge invariance, which is not fixed by the Landau gauge condition (14). However, these 
zero modes occur also for the "unperturbed" operator M = —d 2 and should hence cancel in the ratio Det Mj Det M , 
which we will consider here. 

To illustrate the troubles of the Gribov region in the continuum theory, consider constant gauge fields in SU{2). 
For constant gauge fields the eigenvalues of the Faddeev-Popov operator are given by (see Ref. [11]) 



- 6 Ikl 



a = 0,±l 



(17) 



where 



b=^{Afh){Afk ) 



k = k/|k| 



(18) 



It is seen that even for very small A1 = const there exists always momenta |k| < b so that the eigenvalue of Eq. (17) 
with a = 1 becomes negative and the constant field configuration is outside the first Gribov region. This shows 
that in the continuum theory even arbitrary small (e.g. constant) gauge fields are outside the first Gribov region. 
On the lattice, constant gauge fields are inside the first Gribov region as long as b (18) is smaller than the smallest 
non-zero momentum (the k = eigenvalue cancels against the k = eigenvalue of Mq — —d 2 ). To guarantee that 
the sufficiently small gauge fields are inside the first Gribov region we will introduce an IR-cutoff \i into the proper 
time integral by inserting the factor exp(— ^i 2 r) into Eq. (13). 

Choosing the unperturbed part of the Faddeev-Popov operator as 



M = -d 2 

the free heat-kernel (8) in d-Euclidean dimensions is given by 



(19) 



<z|tfo(T)|i/> = (47rT) 



-d/2 



e 4t 



(20) 



Once the unperturbed part Mo is specified, a recursion relation for the heat coefficients is derived from the heat 
equations (5), (9). With the above choice of M (15) and M (19), inserting the ansatz (7) into the heat equation (5) 
and exploiting the unperturbed heat equation (9) one obtains the following differential equation for the interaction 
part 



d T - d 2 + -(x - y)d + ^-(x - y)A - Ad 
t It 



(x\H(r)\y) =0. 



Using the Landau gauge condition (14) this equation can be expressed as 

d T + -{x-y)D-^-(x-y)A-Dd (x\ H(t) \y) = . 

T AT 



(21) 



(22) 



Inserting here the expansion of the interaction kernel in terms of the proper time (11) one derives the following 
recursion relation for the heat coefficients h n (x,y) 



(/■ + U + l.'- - U)(D-±A 



h k +i(x,y) + {-Ddh k (x,y)) = 



(23) 



with the initial condition 



(x-y)[D--A)h o (x,y) = 0. 



(24) 
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Note that the equations (23), (24) are independent of the number of dimensions so that the heat coefficients h k are 
the same in all dimensions. The dependence on the number of dimensions is entirely contained in the free heat-kernel 
(20). 

Due to the boundary condition (12) the term with n = yields the determinant of the unperturbed kernel (19). 
We therefore find with the explicit expression for the unperturbed kernel (20) from (13) 



(25) 



where 



Trh n = J d d x trh n (x,x) . (26) 
From the representation (25) it is seen that the terms with 

« - \ < (27) 

are ultraviolet divergent, i.e. for n < | we obtain power divergences A d ~ 2 ™ while for n = ^ we obtain a logarithmic 
UV-divergence. Furthermore for fi = the terms with n — | > are infrared divergent, which is a manifestation of 
the presence of gauge fields escaping the first Gribov region as discussed above. Keeping \x finite yields 

where 

T(a,z) = / dTT^e" 1 (29) 

J z 

is the incomplete Gamma function. In the UV- finite terms n > | we can let A — > oo. Then the /i-dependence of these 
terms is given by 

fi d ~ 2n , (30) 

i.e. these terms are diverging for fi — > 0. We are interested here only in the UV-divcrgcnt part of the Faddeev-Popov 
determinant in d — 3, 4, which can be easily extracted. For this purpose it is sufficient to calculate the heat coefficients 
up to including h 2 (x,x), which we will do next. 



IV. CALCULATION OF THE HEAT COEFFICIENTS 



The heat-kernels of generalized Laplacian operators were studied by a number of authors, see e.g. Ref. [12]. Below 
we explicitly calculate the leading and next-to-leading order heat coefficients of the Faddeev-Popov operator (15). To 
this purpose, it is useful to introduce 

D = d+ l -A = D-\A, (31) 

which is essentially a covariant derivative with a rescaled field. The Faddeev-Popov operator can then be re-expressed 
as 

-Dd=-{D+\A)(D-\A) = -D 2 + \A 2 , (32) 

where we have the Landau gauge condition OA = [D, A] = 0, and the recursion relation for the heat coefficients 
Eq. (23) and the initial condition Eq. (24) become 

(x-y)Dh Q (x,y)=0, (33a) 
[(k + l) + (x-y)D]h k+1 (x,y)= [D 2 - \A 2 ]h k {x,y). (33b) 
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Putting k — in the recursion relation (33b) and taking the coincidence limit y — > x yields 

h 1 (x,x)= [D 2 h (x,y)] x=y -±A 2 . (34) 

Acting with D^D V on the initial condition Eq. (33a) yields 

= (D^D U + D u D^)h (x 7 y) + 0{x - y) 
= {2D^ + [D v , D„\ }h {x, y) + 0{x - y) (35) 

In the coincidence limit this yields 

[D li D v h (x,y)] x=y = ±F llv , (36) 

where F M „ is the field stregth tensor constructed with the covariant derivative D (i.e. with the rescaled gauge field 
A/2). This implies that the first term in Eq. (34) vanishes, and the first heat coefficient reads 

tn(x,x) = ~ A 2 . (37) 

With the generators T a in the adjoint representation satisfying trT a T b — —Nq 5 ab we eventually obtain 

Trhx = ^ j ' A d xAl{x)Al{x). (38) 

To find the second heat coefficient h 2 we put k = 1 in the recusion relation (33b), yielding in the coincidence limit 

1 



2 h 2 (x, x) 



D 2 - -A 2 

4 



h 1 {x,y)\ x=v ^(D 2 h 1 {x,y)) x=v + ^{A 2 ) 2 . (39) 



where we have used Eq. (37). To calculate the right-hand side of Eq. (39) we act with the operator D 2 on the recursion 
relation (33b) with k = 

D 2 [l + (x- y)D] h^x, y) = [3D 2 + 0(x - y)] h^x, y) ^ D 2 [D 2 - \A 2 } h (x, y) (40) 
which yields from Eq. (39) 

h 2 (x,x) = l(D 2 D 2 - 1 -D 2 A 2 )h (x, y )\ x=y + ^(A 2 ) 2 . (41) 

The second term in Eq. (41) can be rewritten as 

D 2 A 2 = [Dp, [Dp, A 2 ]] +2[D^A 2 ]D, + A 2 D 2 . (42) 

The last two terms give a vanishing contribution when acting on ho in the coincidence limit (see above), while the 
first one, being a total commutator, gives no contribution when the trace in color space is taken. Thus the whole 
term can be discarded. Applying the operator D 2 D 2 on the initial condition (33a) one gets 

= D 2 D 2 [(x - y)Dh (x, y)] = 2(D 2 D 2 + D fJ _D 2 D fJ _)h (x, y) + 0(x - y). (43) 

Using the identity 

D 2 D 2 - D^D 2 D tl = F^F^ + F M „] D v (44) 
applied to h and adding this to Eq. (43) we get finally 

2D 2 D 2 h (x,y) = (F 2 V + [D^F^DJhofay). (45) 

The commutator term in Eq. (45) gives no contribution in the coincidence limit, since D^hoix, y) = 0(x — y) by the 
initial condition (33a). The second heat coefficient reads then 

\xhi{x,x) = ^tri%J%, + ti{A 2 ) 2 . (46) 
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The first term in Eq. (46) can of course be expressed in terms of the actual field strength tensor F^ v , which is related 
to by 

F IU , = ±F IU ,-±[A II ,A V \. (47) 

The heat coefficient (46) then becomes 




THE COUNTERTERMS FOR d = 3 + 1 COULOMB GAUGE AND d = 4 LANDAU GAUGE 



We arc interested here in the divergent part of the Faddeev-Popov determinant in d = 
start with the d = 3 dimensional case, which corresponds to Coulomb gauge in d = 3 
us. In this case the only UV-divergent term in (25) is n = 1. Since 



= 2 



- r 



this term is IR-finite so that we can take the limit fi 

Det(-Dd) 



In 



Det(-d 2 ) 



div 



yielding 
4tt J 



1 



3 and d = 4 dimensions. We 
- 1 and is of most interest to 



(49) 



3/2 



2A Tr hi . 



(50) 



Inserting the explicit expression for the heat-coefficient hi (38), we obtain for the divergent part of the Faddeev-Popov 
determinant 



In 



Bet(-Dd) 



Dct(-<9 2 ) 



div 



x(A?(x)Y 



(51) 



This result is consistent with Ref. [8] , where it was found that (to the order considered) the renormalization of the 
Faddeev-Popov determinant requires a countcrtcrm of the type (2). Here we have proven that, in general, in d = 3 
there is only one UV-divergent counterterm required and that the coefficient of this term is linearly divergent. 
In d = 4 the terms in (28) with n = 1, 2 are UV-divergent yielding 



In 



Bet(-Dd) 



Dct(-<9 2 ) 



1 



div 



(4tt) 2 



A 2 

A 2 Tr hi + In Tr h 2 
A 1 



(52) 



The heat coefficients h n are the same as in d = 3 and are given by Eqs. (37) and (48). 

A comment is here in order: It is well known that gauge invariance forbids the occurrence of quadratically divergent 
mass terms. In fact, these terms disappear due to a cancellation between contributions from ghost and gluon loops. 
Since we are dealing here only with the Faddeev-Popov determinant, i.e. with the ghost loop, such terms will occur. 



VI. SUMMARY AND CONCLUSIONS 



We have performed a heat-kernel expansion of the Faddeev-Popov determinant in Landau/Coulomb gauge. We 
have shown that in d = 3 Landau (d = 3 + 1 Coulomb) gauge there is a single, unique UV-divergent counterterm 
required, which is given by Eq. (51). It represents a mass term whose coefficient is linearly UV-divergent. In d = 4 
Landau gauge there are two UV counterterms required, which are respectively quadratically and logarithmically UV- 
divergent. The latter term and all the UV-finite terms are IR-singular, reflecting the absence of a mass scale in the 
original theory. However, as is well known the scale anomaly gives rise to dimensional transmutation which manifests 
itself in gluon condensation and the generation of a dynamical mass scale. To access the UV-finite terms in the 
heat-kernel expansion, this dynamical mass scale has to be included in the unperturbed heat-kernel. In principle, 
the heat-kernel expansion is a type of gradient expansion, which should be applicable in the infrared, provided the 
unperturbed kernel is adequately choosen. This expansion is, however, bound to fail near the Gribov horizon, where 
the Faddecv Popov operator develops a zero eigenvalue. 
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